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Abstract 



A systematic investigation of Symanzik improvement in the gauge field action is per- 
formed for the static quark potential in quenched QCD. We consider Symanzik improved 
gauge field configurations on a 16 3 x 32 lattice with a relatively coarse lattice spacing of 
0.165(2) fm. A matched set of Standard Wilson gauge configurations is prepared at (3 = 5.74 
with the same physical volume and lattice spacing and is studied for comparison. We find 
that, despite the coarse lattice spacing, the unimproved and less-expensive Wilson action 
does as well as the Symanzik action in allowing us to extract the static quark potential 
at large qq separations. We have considered novel methods for stepping off-axis in the 
static quark potential which provides new insights into the extent to which the ground state 
potential dominates the Wilson loop correlation function. 
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I. INTRODUCTION 



While string breaking in the static quark potential of full dynamical-fermion QCD has 
long been predicted, direct observation of string breaking in lattice QCD is proving to be 
elusive [|l|-|7| • The main reason for this appears to be due to difficulty in isolating the ground 
state of the static quark potential at large qq separations J7|. Overlap of standard operators 
for separating a static qq pair with excited states of the potential and close spacing in the 
spectrum of the potential at large separations demand significant Euclidean time evolution 
in order to isolate the ground state [|7|-||. 

Accessing large Euclidean times is extremely difficult. APE smearing ]IU| is widely 
recognized as an effective way to approach this region. The smearing of the spatial links of 
the lattice mocks up the flux tube joining two static quarks in the ground state potential. 
The smeared operator provides better overlap between the vacuum and the ground state 
potential and improves the signal to noise ratio in the correlation function [TT|. 



As one approaches the large Euclidean space-time regime, statistical errors grow expo- 
nentially. While it is easy to fit the lattice data for the effective potential to a plateau 
ansatz, it is difficult to have confidence that the asymptotic value has been reached. With 
exponentially growing error bars, correlated x 2 an d goodness of fit parameters can provide 
a lower bound on the Euclidean time regime, but do not provide information on whether 
such a bound is sufficient to isolate the ground state potential. 

Techniques for evaluating the extent to which the ground state dominates the Wilson 
loop are needed. Fortunately for the standard Wilson gluon action, methods exist. However, 
these methods break down when improved actions are used. 

Recently, novel ideas have been explored in the search for string breaking [[7,12-rj]. We 
are motivated by the encouraging results of Ref. [|7| studying string breaking via Wilson 
loops in 2 + 1 dimensional QCD. There, using improved actions on coarse lattices in three 
dimensions, string breaking was observed. These authors emphasize that Euclidean time 
evolution of the order of 1 fm is required to isolate the ground state potential. The efficiency 
afforded by the use of coarse lattices is argued to be key to achieving this goal. 

Here we explore a systematic comparison of the static quark potential obtained from 
standard Wilson and Symanzik improved field configurations. Given the same number of 
configurations, similar lattice spacing and equivalent analysis techniques, we illustrate how 
the Standard Wilson action does as well as the Symanzik improved action in extracting 
the long-range qq potential. While the authors of Ref. emphasize the efficiency of the 
improved action approach, we find that unimproved actions on reasonably coarse lattices 
offer the most efficient and suitably accurate use of limited computer resources. 

Section |TT] briefly explains our technique for calculating the static quark potential from 
Wilson loops. Here we present an alternative way of exploring the off-axis static quark 
potential. This new method provides additional information on the extent to which the 
ground-state potential dominates the Wilson loop and can provide confidence that sufficient 
evolution in Euclidean time has occurred. Section [TTT] outlines the Symanzik improved action. 
In Section [TV], we present and discuss our results and conclude in Section [V|. 
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II. THE STATIC QUARK POTENTIAL 



The spectrum of the static quark potential is determined from Wilson loops W(r,t) of 
area r x t, 

W(r,t) = J2Q(r) exp(-Vi(r)t). (1) 

i 

In order to enhance Ci(r), which measures the overlap of the loop with the ground state 
potential, the spatial links are APE smeared [ IDfl . This smearing procedure replaces a spatial 



link Un(x), with a sum of the link and a times its spatial staples: 

a 3 

U^x) -> (1 - a)Up + 7 E [Uu(x)U„(x + va)Ut{x + fia) 

u=l 

+ Ul(x — va)Un(x — ua)U v (x — va + ixa) . (2) 

This is applied to all spatial links on the lattice followed by projection back to SU(3), and 
repeated n times. 

Tuning the smearing parameters is key to the success of the approach. They govern both 
the Euclidean time extent of the Wilson loop correlation function and the relative contribu- 
tions of ground to excited state contributions. It is vital to have a quantitative technique for 
tuning these parameters in order to optimise future studies of string breaking. The measure 
must be the extent to which the ground state dominates the correlation function. 

Efficient methods exist for the unimproved Wilson action for fine tuning the smearing 
parameters to provide optimal overlap with the ground state potential. For t = 0, W(r, t = 
0) = 1 providing the constraint Y^,iCi(r) = 1 for a given r. For unimproved actions, where 
the transfer matrix is positive definite, each Cj(r) > 0. This means Ci(r) can be monitored 
at large r but small t as the number of smearing sweeps are varied, with the optimal amount 
of smearing occurring when C\{r) ~ 1. The proximity of C*i(r) to 1 for small t may be easily 
estimated from the ratio 

W t+1 (r,t)/W t (r,t + 1) (3) 
which equals C*i(r) in the limit C*i(r) — > 1. This provides a quantitative measure of ground- 



state-dominance for unimproved Wilson actions. We note that it is sufficient [15| to fix the 
smearing fraction, a, and explore the parameter space via the number of smearing sweeps, 
n. 

This procedure can be repeated for a number of alternate paths of links for a given 
separation r. By using variational techniques as described in Ref. |TBJ, the combination of 



paths that gives the greatest overlap with the ground state can be found. 

Actions improved in the time direction do not satisfy Osterwalder-Schrader |T7| positivity 
T8fl . This spoils the positive definite nature of the transfer matrix and the constraint C{(r) > 



is lost. Hence one needs either a new quantitative measure for evaluating ground-state- 
dominance, or to check the merits for using Symanzik improved gluon actions for the static 
quark potential at large r. 



3 



For exploratory purposes, we fixed a = 0.7 and considered the values n = 10, 20, 30 
and 40. The best results on a 16 3 x 32 lattice with lattice spacing pa 0.17 fm, are obtained 
using a = 0.7 and n = 10 for the unimproved Wilson gauge action. This corresponds to a 
transverse RMS radius for the smeared links of 0.31 fm, where the transverse RMS radius 
after n sweeps is defined by 



E,x 2 v; 2 (x) 

E, V 2 {x) 



(4) 



where 



and 



Vi{x') = £ 



[1 - a)5 x t x + — ^2{S x > tX+fl + 5 X > )X . 

4 M=l 



V (x) 



K-i(^) 




Analogous to the RMS smearing radius of Jacobi Fermion Source Smearing ]T{| , this provides 
a reasonable estimate of the transverse smearing of mean field improved links (U^(x) ~ 1). 

Wilson loops, W(r, t), or more precisely W(x, y, z, t) where r 2 = x 2 + y 2 + z 2 , are calcu- 
lated both on-axis, along the Cartesian directions, and off-axis. On-axis Wilson loops are 
those that lie, e.g., in the x — t plane only; off-axis loops begin, for example, by first stepping 
into the y or z (or both) directions before proceeding through the x — t plane.. This provides 
an alternative to the usual method of calculating the off axis potential by building paths in 
three different directions using small elemental squares, rectangles or cubes and multiplying 
them together to form larger paths (see, for example, Ref. @]). These standard techniques 



for calculating the off-axis potential may be combined with the approach described here 
using the variational method described extensively in Ref. . 



Due to the periodicity of the lattice, the size of our Wilson loops are limited from 1 to 
a little over half the smallest lattice dimension in the on-axis directions and between and 
3 in the transverse directions. For example, for a 16 3 x 32 lattice, the sizes of the Wilson 
Loops, t x x x y x z, vary from alxlxOxO loop, to a 10 x 10 x 3 x 3 loop. Statistics 
are improved by transposing the loops over all points on the lattice and by rotating through 
the three spatial directions 

In order to efficiently calculate Wilson loops of various sizes, including off-axis loops, we 
build products of links in each direction that we are considering for our Wilson loop. Link 
products extending from every lattice site are calculated in parallel. The loops are then 
formed in parallel by multiplying the appropriate sides together. Each side is created by 
reusing the components of the previous loop and one additional link. The same approach 
can be extended to loops that travel off-axis. 

In an attempt to isolate the ground state potential for off-axis paths, one symmetrizes 
over the path of links connecting the off-axis heavy quark propagators, exploiting the full 
cubic symmetry of the lattice. For example, two points separated by n x sites in the IE- 
direction, n y sites in the ^-direction and n z sites in the ^-direction, may be connected by 
n x link products in the x-direction, n y link products in the ^-direction and n z link products 
in the ^-direction which we denote by the triplet xyz. Instead of only calculating off- 
axis paths in the specific order xyz, we average over spatial paths calculated in the order 
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xyz, xzy, yxz, yzx, zxy, zyx. We calculated loops using this path-symmetrized technique as 
well as loops using a non-path-symmetrized operator where the order xyz alone is considered. 
The former form of operator is designed to suppress excited states by incorporating the full 
hyper-cubic symmetry of the lattice, whereas the latter operator is susceptible to excited 
state contamination. By comparing the static quark potential for these two operators, one 
can gain qualitative information on the effect of excited states in the static quark potential. 

III. NUMERICAL SIMULATIONS 



The tree-level 0(a )-Symanzik-improved action |2T| is defined as, 



So = ^fE^Tr(l - U sq (x)) - -^5>eTr(l - U rcct (x)), (5) 

sq llU rect 

where the operators U sq (x) and U rect (x) are defined as follows, 

U sq (x) = U^x)U v {x + p)Ul(x + u)Ul(x) (6) 
U Iect (x) = U^{x)U v {x + fi)U v {x + v + fi)Ul(x + 2v)Ul(x + O)Ul(x) 

+ U^U^x + jl)U v {x + 2fi)Ul(x + /t+ 9)Ul(x + u)Ul(x). (7) 

The link product U rcct (x) denotes the rectangular 1x2 and 2x1 plaquettes. u is the 
tadpole improvement factor that largely corrects for the quantum renormalisation of the 
operators. We employ the plaquette measure 

n \ 1/4 

u = (-neTr(U sq )) . (8) 

Eq.(^|) reproduces the continuum action as a — > 0, provided that (3 takes the standard value 
of 6/ ' g 2 . Perturbative corrections to this action are estimated to be of the order of two to 
three percent p2 



Gauge configurations are generated using the Cabibbo-Marinari f23| pseudoheat-bath 



algorithm with three diagonal SU(2) subgroups. Simulations are performed using a parallel 
algorithm with appropriate link partitioning p4j . Configurations are generated on a 16 3 x 32 



lattice at (3 = 5.74 using a standard Wilson action which corresponds to lattice a spacing 
a = 0.165(2), and on a 16 3 x 32 lattice at (3 = 4.38 using the Symanzik improved action 
(Eq. 0) which also corresponds to a lattice spacing a = 0.165(2) fm. Thus the two lattices 
have the same lattice spacing and physical volume. 

Configurations are selected after 5000 thermalization sweeps from a cold start. The 
mean link, uq, is averaged every 10 sweeps and updated during thermalization. For both 
the standard Wilson action and the Symanzik improved action, configurations are selected 
every 500 sweeps. The following analysis is based on an ensemble of 100 configurations for 
each action. 
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IV. SIMULATION RESULTS 



The effective potential is obtained from 

which is expected to be independent of t for t ^> 0. Figure [L] displays the effective potential as 
a function of Euclidean time, t, obtained from 100 configurations generated via the Symanzik 
improved action at (3 — 4.38. For r > 7, we find that the signal is generally dominated by 
noise for t > 4, so we set the upper limit of our fitting range to t max = 4. The good plateau 
behavior at small Euclidean time is a reflection of the optimized smearing. Choosing the 
lower limit t m i n = 1 leads to large x 2 /d.o.f for large r. We fix the fitting range to be, in most 
cases, t = 2 to 4. The string tension is then extracted from the ansatz, 

V(r) =V + ar- e/r (10) 

where e = n/12 pB| , and V and a are fit parameters. 

The values for V(r) + e/r = ar + V are then fitted to a line and the slope, a, and 
intercept, V , are extracted. We fit the range 3 < r < 7 such that we are not sensitive to 
the coulomb term and its associated dicretisation artifacts |26|| . The error analysis is done 
using the Bootstrap method and all errors quoted are statistical only. Using y/a = 440 MeV 
to set the scale, the lattice spacings are determined. 



A. Wilson Loop Correlation Function 

To have any hope of seeing string breaking in QCD with dynamical fermions, the gauge 
links must first be smeared to improve the overlap of the ground state. The effects of 
smearing the spatial links (Eq. (|2])) before calculating Wilson Loops is well known. Smearing 
provides access to the static quark potential at larger distances, providing a better chance 
of eventually being able to detect string breaking. Also, the lines with smeared links exhibit 
better plateau behavior in V t (r) than the unsmeared ones, indicating that we have better 
isolation of the ground state. 

Figs. [I] and |^ illustrate the time dependence of Eq. ^|for Symanzik improved and standard 
Wilson gluon actions respectively. For clarity, we only plot values for r that are obtained 
from on-axis loops. For both the improved and unimproved actions, we see a clean signal 
up to r « 8 for the first three or four time slices. For r values grater than this, the quality 
of the statistical signal does not allow values of t larger than 2 or 3 to be included in the 
fits. 



B. Off Axis Perturbations and Symmetry 

Our method for stepping off-axis requires the use of a path-symmetrized operator. Fig. |3| 
displays results for the path-symmetrized operator for separating the qq pair, while Fig. |] 
illustrates the non-path-symmetrized result. Both plots are obtained with 10 smearing 
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FIG. 1. Euclidean time evolution of the effective potential. From bottom up, the "horizontal" 
lines correspond to r = 1 through 10. Here we use the Symanzik improved action with (3 = 4.38 
and 10 sweeps of smearing at a = 0.7. Axes are in lattice units. Points are offset to the right for 
clarity. 
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FIG. 2. Euclidean time evolution of the effective potential. From bottom up, the "horizontal" 
lines correspond to r = 1 through 10. Here we use the standard unimproved Wilson action with 
(5 = 5.74 and 10 sweeps of smearing at a = 0.7. Axes are in lattice units. Points are offset to the 
right for clarity. 
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FIG. 3. The static quark potential, V(r), as a function of the separation r. Data is from 
the Symanzik improved action at /3 = 4.38 with 10 sweeps of smearing at a = 0.7 and a 
path-symmetrized operator. Time slices t = 1 to 4 are used in the fit of the correlation func- 
tion. Error bars are magnified by a factor of 200 and placed on the r-axis for clarity. 



sweeps at a = 0.7. Banding is clearly evident in the non-path-symmetrized case but as soon 
as we path-symmetrize our operators, the banding largely disappears in off-axis points up to 
about r = 7. After r = 7, the banding still persists in off-axis points where loops of the form 
rxlx3,rx2x3orrx3x3 (where r = on-axis step size and is greater than 7) are used. 
However the banding is negligible in off-axis points up to, and including, r x 2 x 2. Since 
our off-axis points are obtained by firstly stepping in one direction and then stepping in a 
different Cartesian direction, we are forming a right angle, not an approximate straight line. 
The fact that for r < 7, the r x 3 x 3 points lie on the line-of-best-fit is an interesting result 
and provides support for the use of this operator in a variational approach when searching 
for ground state dominance and its associated string breaking in full QCD. However, the 
persistence of banding at large r is an indicator that further Euclidean time evolution is 
required to isolate the ground state potential. 



C. Effect Of Using Symanzik Improved Gauge Field Configurations 

Historically, the main feature of improved gauge field actions is the improved rotational 
symmetry ||22|| . This can be seen in our configurations by comparing Figs. [5] and |6|. These 
graphs are enlargements of the small-r area. The off axis points for the Symanzik improved 
P = 4.38 lattices lie closer to the line of best fit through the Cartesian points, r = 3 to 7, 
than the unimproved Wilson (3 = 5.74 lattices. 

However, here we are most interested in the large distance properties of the Wilson loop. 
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FIG. 4. The static quark potential, V(r), as a function of the separation r. Data is from 
the Symanzik improved action at /3 = 4.38 with 10 sweeps of smearing at a = 0.7 and a 
non-path-symmetrized operator. Time slices t = 1 to 4 are used in the fit of the correlation 
function. Statistical error bars are too small to be seen. 



To draw firm conclusions on the merits of using Symanzik improved actions in the search 
for string breaking, we create a matched set of standard Wilson gluon configurations, tuned 
to reproduce the string tension of the (3 = 4.38 improved configurations. 

Figs. | and compare Symanzik improved and unimproved Wilson results for V(r) 
extracted from correlated fits within the range 1 < t < 4. A similar comparison is made 
with Figs. § and |9| within the range 2 < t < 4. Magnified error bars are plotted on the X-axis 
to allow a comparison of signal to noise between the two actions. Exponential growth in the 
error bars for large r is apparent in the fits with 2 < t < 4. 

It is crucial to compare matched lattices with the same lattice spacing and same physical 
volume. For example, a comparison of an unimproved Wilson lattice with (3 = 5.70 and 
lattice spacing a = 0.181(2), with a Symanzik improved lattice at (3 — 4.38 and a = 0.165(2) 
leads one to incorrectly conclude that Symanzik improved actions lead to an improvement in 
the signal-to-noise ratio of the static quark potential at large qq separations. The unimproved 
configurations with (3 = 5.70 lose the signal around r/a = 7.5 (r = 1.37 fm) when using 
time slices t = 2 — 4, whereas the improved configurations hold the signal up to r/a = 9.5 
(r = 1.59 fm) when using t — 2 — 4. This effect is due to the slightly larger lattice spacing in 
the (3 = 5.70 simulations, which spoils the signal-to- noise ratio even after the larger lattice 
spacing is taken into account. 

Close inspection of Figs. [7| and ||, where the Euclidean time regimes 1 < t < 4 and 
2 < t < 4 are compared, reveals how the off-axis points can also be used to gain confidence 
in ground state dominance. In Fig. 0, banding is apparent. However, upon adjusting the 
Euclidean time regime for the fit to larger times in Fig. || the banding is largely removed. 
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FIG. 5. A close up of the static quark potential, V(r), as a function of the separation r. Data 
is from the Symanzik improved action at f3 = 4.38 with 10 sweeps of smearing at a = 0.7 and a 
path-symmetrized operator. Time slices t = 2 to 4 are used in the fit of the correlation function. 
The dashed line is a fit of Eq. |l(] to on axis points r = 3 to 7. 
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FIG. 6. A close up of the static quark potential, V(r), as a function of the separation r. Data 
is from the unimproved Wilson action at /3 = 5.74 with 10 sweeps of smearing at a = 0.7 and a 
path-symmetrized operator. Time slices t = 2 to 4 are used in the fit of the correlation function. 
The dashed line is a fit of Eq. pi] to on axis points r = 3 to 7. 
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FIG. 7. The static quark potential, V(r), as a function of the separation r. Data is from 
the unimproved Wilson action at j3 = 5.74 with 10 sweeps of smearing at a = 0.7 and a 
path-symmetrized operator. Time slices t = 1 to 4 are used in the fit of the correlation func- 
tion. Error bars are magnified by a factor of 200 and placed on the r-axis for clarity. 



V. CONCLUSIONS 

We have calculated the static quark potential in quenched QCD using Symanzik im- 
proved and unimproved Wilson gluon actions. We have kept the lattice spacing and the 
physical volume of these lattices equal so that we can meaningfully compare the results. 
The number of gauge field configurations (100 here) is also held fixed for each action. We 
have explicitly shown that, despite the relatively coarse lattice spacing, the unimproved and 
computationally less expensive Wilson action does just as well as the improved action in ex- 
tracting the qq potential at large separations. If one wishes to keep non-perturbative physics 
such as non-trivial topological fluctuations on the lattice, then one needs a < 0.15 fm p7| . 



and thus r/a > 7. In this case, the unimproved, standard Wilson gauge action is ideal for 
todays string breaking searches and computational resources can be redirected elsewhere. 
Another advantage for using unimproved Wilson gauge configurations is that we recover the 
extremely useful method for calculating the overlap with the ground state, C\(r), and thus 
tuning the smearing parameters. 

We also explored the use of unconventional paths in accessing off-axis values of r in the 
static quark potential. These paths can provide insight into the extent to which the ground 
state potential dominates the Wilson loop at large Euclidean times. Provided the paths are 
symmetrized, these new paths provide useful information on the ground state potential and 
nearby excited potentials. Combined with standard paths and variational techniques, these 
paths offer additional promise for the search for string breaking in lattice QCD. 
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FIG. 8. The static quark potential, V(r), as a function of the separation r. Data is from 
the Symanzik improved action at (3 = 4.38 with 10 sweeps of smearing at a = 0.7 and a 
path-symmetrized operator. Here time slices t = 2 to 4 are used in the fit of the correlation 
function. Error bars are magnified by a factor of 200 and placed on the r-axis for clarity. 
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FIG. 9. The static quark potential, V(r), as a function of the separation r. Data is from 
the unimproved Wilson action at j3 = 5.74 with 10 sweeps of smearing at a = 0.7 and a 
path-symmetrized operator. Time slices t = 2 to 4 are used in the fit of the correlation func- 
tion. Error bars are magnified by a factor of 200 and placed on the r-axis for clarity. 
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